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INTRODUCTION 



The problem of ground effect has "been discussed from 
time to time in this periodical, the most recent instance 
"being in an article "by the author (reference 24) wherein 
he analyzes the problem on the basis of his developed bi- 
plane theory. The conclusion of that article is that the 
increment of lift due to ground effect is largely attrib- 
utable to the effect of induction of the free vortices, 
and is practically equivalent to a virtual increase in 
aspect ratio. The ground clearance was of the order of 
magnitude comparable to the wing chord. 

Recently we have received some new articles and stud- 
ies on this interesting subject: an inf ormat ive^ar t icle 
by M. Le Sueur (reference 31), a report by G. Datwyler 
(reference 22), which treats the case of minimum distance 
from the ground and is confined to the plane problem only, 
and lastly, some theoretical studies by Tomotika and oth- 
ers (references 23 and 25), also confined to the plane 
problem. 

We shall briefly review these reports without regard 
to chronological order. 



I. LE SUEUR ! S REPORT 



Ground effect is a very controversial subject, both as 
to cause and to its effect. All observations made in free 
flight are in agreement as regards the existence of cer- 
tain systematic phenomena such as the greater facility of 
low-wing airplanes to take off, impossibility of certain 



* M Teorie ed esperienze sul problema dell ? ala in vicinanza 
del suolo. From Pubblicazioni della R. scuola d f In- 
gegneria di Pisa, series VI, no. 261, July 1935, pp. 
1-25. 
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heavily loaded long-distance airplanes to gain altitude, 
the prolonged glide of low-wing airplanes on landing, etc* 
Even so, there exist erroneous and misleading opinions as 
to the true cause of these phenomena. 

It is a certainty that ground interference lowers the 
drag, assuming the lift to "be equal, and in quite notice- 
able proportions. As to the maximum lift, there is no 
theory to attest to its increase; in fact, divers experi- 
ments in accordwith certain theories appear to indicate 
occasionally, a decrease. 

On these premises we shall examine the various phases 
of the problem from the theoretical as well as from the 
experimental point of view. 

A. THEORIES ON INTERFERENCE EFFECT 



These theories are able to treat the plane problem 
(wing of infinite span) or even the three-dimensional 
problem. In any case, the study on interference reverts 
to the general study of the biplane making use of the 
artifice known as "reflection. ,r There are available for 
this purpose (three-dimensional problem) the noted Prandtl- 
Betz formulas, which may be applied in our case and which 
were published by Wi esel sberger (reference 5). These for- 
mulas - " were chosen in 1924 by Toussaint for comparison with 
his experimental values of the coefficient in different 
cases of monoplanes and biplanes with ground interference. 
They are generally accepted as sufficiently approximate, 
though studies made subsequently disclose results more or 
less at variance with the former. 

A report by Rosenhead (reference 18) treats the lift 
of a flat plate between parallel walls - an analysis based 
on a method of conformal^ transformation , and which he com- 
pares with Glauert's results (reference 32). 

There is also a study by Muller (reference 19) in 
which he applies to two symmetrical airfoils Ferrari's 
method of conformal transformation, the results of which 
are at variance with experiment, as they are indicative of 
a decrease in lift. 

Finally, there is the study by Pistolesi (reference 
24) who, in the case of the wing with infinite span, finds 
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that the circulation increases with the angle of incidence 
up to a certain value of this incidence, "beyond which it 
decreases again; and, taking into consideration the change 
in horizontal speed he expresses the change in lift coef- 
ficient due to interference "by simple formulas (reference 
24). He then passes to the case of finite span and finds 
the percent change of lift in function of the incidence 
for divers values of h/l (I = chord) - which results are, 
at least qualitatively, in agreement with experiment. In 
particular, the increase in lift decreases until it is 
canceled out as the incidence increases. 

Lastly, there are the Japanese reports, the final 
conclusions of which are identical with those of Pistolesi: 
At low incidence the lift increases as the distance from 
the ground decreases; at high incidence, on the other hand, 
the lift decreases as the distance from the ground de- 
creases. 

B. EXPERIMENTAL METHODS 



aJ_Teji^_wi.ih — .smal lrscale_models,- Not wishing to go 
"back as far as Betz ! s experiments of 1912 (reference l) 
(which, while disclosing negligible interference effects, 
were quite inaccurate), we record an interesting study "by 
Cowley and Lock (reference 3), "based on tests made in Eng- 
land, in July 1920, in the 4-foot No. 1 wind tunnel, at 
40 feet per second wind speed, on an R.A.F. 15 "biplane of 
zero stagger, in connection with the "Tarrant" "biplane. 

The ground was represented, in part ("by real ground) 
by a vertical sheet of tin 4 feet high, 3 feet long; and in 
part by two identical biplanes symmetrically disposed, in 
respect to virtual ground, according to the reflection 
method. Measurements were made of the lift, drag, and 
pitching moment at angles of attack ranging from -6° to 
14°, and for values of h/b (b s span) ■ 0.167 and 0.306. 

At about the same time, the Massachusetts Institute of 
Technology also made some similar tests in its 4-foot tunnel, 
at wind speeds of 30 miles per hour, except in two cases 
where it was raised to 40 and 45 miles per hour. These 
tests, reported by A. E. Raymond (reference 6), were made 
on three 3- by 18-inch models: a Martin No. 2, an-. R . A f F • 15 
special, and a U.S.A. 27. Experiments were also made by 
the flat-plate method (3*#ly birch 3/8-inch thick, 4 feet 
high, 3 feet wide, with leading edge chamfered on the side 
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away from the model) and "by the reflection method. In 
"both cases the tests were run for different angles of in- 
cidence and different ground distances, varying from l/4 
I to 21 (I = chord). 

The same experimental method was used in 1921 in Ger- 
many to check Wiesel sherger * s formula (reference 5). The 
model was a monoplane of 124 cm (48.32 in.) span, aspect 
ratio 9. 

Some years later Toussaint undertook a series of sys- 
tematic experiments in the 6^-foot No, 1 wind tunnel at 
Saint Cyr (reference 9); the ground was represented hy a 
sheet of aluminum 4 mm (0.157 in.) thick, and 1.60 m (5.24 
ft,) long. The measurements were made on a wire "balance; 
the wind speed was 32 to 33 m/s (105.0 to 108.3 ft. /sec). 
The models were a Liore (S.C. 133 a) wing, a Fokker (S.C. 
106 a) wing, a Fokker (S.C. 106 a+h) biplane wing, as well 
as two Breguet 14 A 2 airplane mo dels of 1/10 and l/ 20 
scale. The tests were made very painstakingly to assure 
accurate results. The incidence was changed in stages of 
3° each from -9° to +15°, and for three ground distances: 
0.530 m (1.74 ft.), 0.438 m (1,44 ft.), and 0.240 m (0.787 
ft.). These test data served as check on Betz's formulas, 
to which we shall refer again later on. 

Other tests were made in the Eiffel wind tunnel on 
models mounted in the presence of a platform representing 
the ground. The wind speed was 25 m/s (82 ft. /sec); the 
model was a Caudron R.220; the distance from the ground 
0, 100, 200, 300, 400, and 500 mm (3.94, 7.87, 11.81, 15.75, 
and 19.69 in.). Unfortunately, the use of only two inci- 
dences, 0° and 12°, detracts from the value of these ex- 
periment s • 

Next, there are the tests of the Wihaul t-Penho et com- 
pany on an airfoil 172, aspect ratio 5, a low-wing mono- 
plane, type 313, airfoil 209, aspect ratio 7.8; a low-wing 
monoplane 280, airfoil 125, effective aspect ratio 6.84, 
real aspect ratio 7.85; a low-wing monoplane 287, airfoil 
215, aspect ratio 8.4; a low-wing seaplane 240, The mod- 
els had spans ranging from 1 meter to 1.30 meters; ground 
distance and incidence were changed systematically. 

Lastly, we mention the tests reported "by Datwyler 

(reference 22), which refer to very short distances from 

the ground, and which had "been conducted in Gottingen 
and Zurich, 
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From among the tests made on aerodynamic carriages, 
those of Tonnies deserve special mention (reference 2l)« 

b_)__Full — scale experiments; ,~ Here we have first of all 
the flights made during 1927 in the United States, as re- 
ported "by Elliott G. Reid, in Technical Report No. 265 
(reference 12). They were made with a Vought VE-7 bi- 
plane whose aerodynamic characteristics including those of 
the propeller, had "been established in previous tests. 
The flights were made with the lower wing at 1.50 to 3.75 
(5 to 9 ft.) above the ground. An interpolation method 
was used to eliminate the necessity of maintaining strict- 
ly level flight. 

Other interesting full-scale tests are those cited in 
Tifnnies 1 report (reference 21), made on a Klemm 26-2a at 
Hanover. The records were made with a motion-picture 
camera at 2, 4, 7, 10, 15, and 20 m (6.56, 13,12, 22.97, 
32.8, 49*2, and 65.6 ft.) height. The tests were numer- 
ous and followed a set schedule. 

The writer of that report offered some suggestions 
regarding the photographic method and advocated a method 
based on the emission of radio signals. 



C. TEST DATA - THEIR AGREEMENT WITH 0I7E ANOTHER 

AND WITH THEORY 



The results of the various experiments can always be 
satisfactorily interpreted with Wi e s el sberger f s formula: 

Aa = - a 

c p 2 s 

AC r , - a — w 

which can be translated in a change in aspect ratio 

A' = 

l - a 



However, there is no accord between the reflection and 

the flat-plate method, especially for very short distances. 
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Figure 1 offers an example of this discrepancy. 

This was noted "by Cowley and Lock in 1921 (reference 
3) who - for the plate - imputed the disturbance set up 
by its leading edge, which causes the air flow to deflect 
upward, while as regards the reflection method, they raise 
the element of doubt about the assumption of symmetrical 
flow which could in reality be replaced by an oscillating 
flow. The authors conclude that ground interference pro- 
duces an appreciable effect on the pitching moment, but 
very little on maximum lift; the maximum L/ D is increased 
from 10 to 13 by the reflection method, and from 10 to 15 
by the flat-plate method. 

For the plate at 38 mm (1.496 in.) height, equivalent to 

Z = T (h = twice the ground distance), the increase in 
b 6 

lift and in L/D far a given angle of attack, is about 
twice as high by the flat-plate method as by the reflec- 
tion method. 

These findings are in agreement with those of other 
experimenters, for instance, Raymond (reference 6). 

The Eiffel tests, made exclusively by the flat-plate 
method disclosed systematically a very much greater influ- 
ence than Wi e s el sberger 1 s formula stipulated, and furnished 
interference factors a which were about twice the theo- 
retical figure. 

Le Sueur 's conclusion is that the wind-tunnel tests 
are in agreement with Wiesel sberger 1 s theory. 

Tonnies 1 carriage experiments (reference 21) are in 
accord with the theory; they record a decrease in lift at 
high incidence as stipulated by the theory (Pistolesi, 
Tomo t ika ) . 

The free-flight tests by Reid (reference 12) are in 
close accord with the theory. This also holds true, in 
general, at least, for the flight tests described by Tun- 
nies (reference 21), 

Le Sueur comes to the conclusion that as yet we >are 
actually unable to make any definite deductions, although 
for great ground distances and small angles, the different 
experiments seem to be in agreement with Wi e sel sberger f s 
theory, which likens the ground interference to a ficti- 
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tious increase in aspect ratio. The effect in flight cor- 
responds to what in the United States is called "float- 
ing. ,r For great angles and short distances from the 
ground there seems to he a loss of lift; perhaps it is the 
phenomenon which is called "pancaking" in the United States* 

Lastly, at very high angles of attack and very short 
distances from the ground, a marked increase in lift may 
result (air-cushioning effect), according to Latwyler f s 
experiments. 

In the next chapter, Le Sueur discusses the effect of 
the phenomenon on the different phases of airplane motion 
close to the ground, comparing high-wing to low-wing air- 
planes from the point of view of ground effect which, how- 
ever, has no direct bearing on the argument which we shall 
advance. Suffice it to say that a lift increase may "be 
useful in any case; the decrease in drag, useful at take- 
off, may become detrimental at landing as it increases 
the landing run, 

II. DATWYLER. r S THESIS ON GROUND EFFECT 



Datwyler starts with some curious experiences. While 
Wiesel sberger 1 s theory leads one to presume a decrease in 
maximum lift (effect of induced horizontal velocity, at 
least in the plane problem), he observes the fact reported 
by Zimmermann in the aer o t echni cal debate at the Polytech- 
nic Institute, Zurich, that in 1932 an American low-wing 
monoplane (Lockheed) landing with retracted landing gear, 
actually had a lower emergency landing speed than expected. 
The pilots reported they felt as if at the last moment be- 
fore touching the ground, an air cushion had formed below 
the wing. The following year, 1933, an airplane of the 
same type actually had to make a forced landing with 
wheels retracted, when the retraction mechanism failed to 
f unc t i on . 

From these occurrences the author concludes that for 
small ground distances the omission of th-e finite profile 
dimension and the substitution of the profile by one sin- 
gle bound vortex of given circulation, is no longer per- 
missible. For wings of finite span the problem becomes, 
of course, more complicated. 

ti 

The theoretical part, detailing the Gottingen and 
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Zurich experiments, is followed by an appendix on the min 
imum induced drag of a wing with dihedral. 

A. THEORETICAL PART 



He "begins with the treatment of the two-dimensional 
flow about a flat plate whose rear edge rests on an infi- 
nite plane (fig. 2) at an angle a. The problem is re- 
solved "by conf ormal transformation of plane z = x + i y 
on the plane t = r + i s (fig. 3), whereby the charac- 
teristic function is a particular case of the Christoffel 
Schwarz differential equation. 



Putting 



t 
c 



= T 



(c may be arbitrarily assumed) and 



b 
c 



with 



0 , 77 e have 



z - z. 



(- t) 



(0 - T ) (1 - t)' 



d t 



Integration along the contour of the plate, with s = di s 
tance of a point of the contour starting from the trail- 
ing edge and I = plate length, we obtain: 



1-p 



in which the author uses the function r* 

He determines the lift of the plate on the following 
considerations: The drag is zero in a perfect fluid. Ac 
cording to figure 4, the lift P resulting from the nor- 
mal pressures on the plate fcp) and the suction force 
S at the plate tip, is equal to S/ sin a. The suction 
force S is substantially conditional to the velocity 
distribution in direct proximity to the plate tip. He 
establishes the lift coefficient at 



,_TT 

sin (ttP) 



(nf-i) 



1-20 



(1) 



as against 



= 2tt sin (ttP) 



(2) 
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for the free plate- Equation (l) for small values of P 
can "be approximated: 

c p = fT 23 

Next the author determines the moment with respect to the 
leading edge at 

M = | U 2 I 2 (J Q - J u ) (3) 
where J Q and J u indicate the integrals 



f CW 0 - f) * (f 

s = o 



- 7 s = o 

and whereby v 0 and v u denote the values of the velocity 

on the top and bottom sides of the flat plate. The values 
of these integrals are plotted. From (3) we deduce: 

C m = J 0 - J u (30 

Another case treated "by this author is the lift of a 
flat plate close to the ground at 90° angle of attack* He 
employs the reflection method and a conformal transforma- 
tion on plane X = $ + £$< The result is expressed in 



P X + 1 \ ay 



while X = — 1 — (4) 

a - "b 

and b and a denote the ordinates of the two edges (upper 
and lower) of the flat plate. 

Lastly, the author proposes an approximate evaluation 
in the general case, "based on the effect of the reflected 
wing. This induces on the reflected wing a horizontal ve- 
locity which the author di|regards as small compared to V, 
and a vertical velocity w in upward direction* He com- 
putes this velocity in the center of gravity of the circu- 
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lation; that is, at l/4 plate chord on the "basis of the 
plotted velocities induced "by the flat plate when moving 
downward at constant speed w. This flow is, up to a small 
share of the circulation, equivalent to the downwash on 
the flat plate. The speed w is that which corresponds 
to the effective angle of attack, or: 



For w* = zw, it is: 



hence 



2 = a + S- 



S = a + *2 

V V 

w = a 

V 1 - z 



and since Cp = 2tt a> e ff B 2tt ? 

= 2tt — &~ (5) 
l J 1 — z 

Thus, we ohtain: 

M =0 0.025 0.05 0.10 0.20 oo 

if 

r L- = oo 6.67 3.51 2.00 1. 389 1 

1*< z 

where s = distance above ground, and I = chord. Such 
result, as the author justly recognizes as valid for small 

angles of attack and therefore — — — , renresents the an- 

1 m Z 1) 
gular coefficient of the tangent of the lift coefficient. 

The author then plots the lift coefficient for various 

s/l (fig. 5). 

B i EXPERIMENTS 



These are also given in two parts. The first experi- 
ments were made at G-ottingen in the small wind tunnel on a 



The last part of this report contains a more exact and 
complete application of this concept. 
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thin, symmetrical airfoil (thickness parameter 6.075) 
evolved from a Joukowski airfoil, as follows: To assure 
a flat "bottom side and, for reasons of symmetry, a flat 
upper side, straights were drawn from the trailing edge 
as tangents to the profile curve, forming a trailing- 
edge angle of exactly 8°. This gave the trailing edge 
the desired strength. The radius of curvature of the 
leading edge was also enlarged. The airfoil had a span of 
800 millimeters, a wing chord of 200 millimeters, and an 
aspect ratio of 4. A flat 20-millimeter thick hoard with 
symmetrically tapered leading edge represented the ground. 
It fitted directly to the nozzle flare and extended with 
a length of 1,560 millimeters, sufficiently far beyond the 
trailing edge of the airfoil. Being 1,300 millimeters 
wide, it extended 50 millimeters at either side of the 
lm 200-millimeter wide air stream. The wind speed was 30 
meters per second; the Reynolds Number 

Re = - 3.7 x 10 5 

The tests were made with and without end plates, and in the 
case of the isolated wing, with proper corrections for au- 
toinduction and scale effect, but disregarded as negligi- 
ble in the case of a wing in proximity of the ground. The 
ground distance, measured at the trailing edge, was 5, 10, 
20, and 40 millimeters. 

Apart from secondary phenomena (negative lift at very 
small angle of attack due to a venturi effect between air- 
foil and plane), the test data were qualitatively in accord 
with the theory: diminution of drag for equal lift; in- 
crease of lift at equal angle of attack. However, the 
lift increase falls below the theoretical expectations. 
Moreover, there is no change when passing from 10- to 5- 
millimeter slot width. The reason for this is to be found 
in the retarding effect of the flat plate used to repre- 
sent the ground which - rather than being immobile with 
respect to the wing - is stationary with respect to the 
fluid. Figures 6 and 7, photographed in the water tank, 
manifest the unlike behavior in both cases and the unlike 
formation of the boundary layer. (Fig, 7 corresponds to 
moving fluid and fixed wing; fig. 8, to moving wing and 
fluid at rest.) To overcome this drawback, the author had 
recourse to a second series of experiments in the small 
wind tunnel of the Federal Technical Institute at Zurich, 
where he used the method of reflection which, although not 
reproducing exactly the real phenomenon (no allowance for 
friction of air entrained by the wing against the ground 

mm x 0.03937 = in. m/s x 3.28083 = ft. /sec. 
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and the cumulative lift which it is necessary to overcome). 
The tests were limited to the determination of" the lift. 
The results are much more in accord with theory than the 
previous ones, as seen on figure 8. 

If, however,- the obtained lift curves are extrapolat- 
ed to zero, it will he noted that the incidence for zero 
lift varies as the distance. This, the author ascribes 
to the effect of finite profile thickness. If the curves 
are shifted approximately 1.8° to the left, so as to afford 
Op = 0 for a = 0, the result will be the curves of fig- 
ure 9. The results by the reflection method are much more 
appropriate than by the flat-plate method. The use of side 
screens tends to prevent the lateral escape of air dammed 
up on the bottom side of the wing and consequently, a lat- 
eral pressure decrease. The lift distribution becomes 
more complete and higher lifts may be expected. 

Pressure measurements in the slot between the two 
wings showed positive pressure by adhering flow and, con- 
sequently, retarded flow. After separation of the flow 
the negative pressure on the upper side of the wing con- 
tinues into the slot itself. This separation occurs soon- 
er as the slot becomes narrower; that is, as the angle of 
attack becomes smaller. 

This interesting report concludes with an appendix, 
in which the author uses the results obtained in the theo- 
retical study to solve the problem of minimum induced drag 
of a wing with dihedral angle (fig. 10). He determines 
the potential, also the potential jump &cp = cp 0 .~ cp u be- 
tween the two sides of the plate by a current having an 
asymptotic velocity equal to unity. This jump multiplied 
by 2w (w = induced velocity (constant)) on the wing sup- 
plies the circulation in the corresponding point of the 
wing, that is, T. There is, according to Kut ta- Joukowski , 
a lift , 




o 



whereby T 0 indicates the circulation at wing center; 
an induced drag 



al so 




r = 2w (cp 0 - cp u ). 
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With \|/ as the ratio "between induced drag R± of 
the wing with dihedral, and that of the wing with ellip- 
tical distribution and equal span (b), we have: 

2 

, cose f „ N 

= 2tt — ( 6 ) 

(it - 2€) (tt + 2c) 

The value \1/ is plotted in figure 11 against e. 

The report closes with a notation on Tomotika's exper- 
iments which pertain, however, to ground distances still 
amounting to 20 percent of the wing chord. 



III. OTHER THEORETICAL STUDIES 
ii 



A) LIULLER (reference 19) 



Muller proposes to effect the conformal transformation 
of two circles symmetrically disposed with respect to axis 
x in two profiles, also symmetrically disposed with re- 
SDect to the same axis. 



His transformation is as follows: 



z = 



+ — 



a 



^ + — "72 



(?) 



The singular points in which 
found within the two circles. 



dz/d£ cancels out must he 
Thus , 



dj2 _ 

H = 



a - £7 



(8) 



Equated to zero, we find the four singular points 

Fx U« 



u ti 



If a is real, we ohtain a "biquadratic equation, and 
that is the case which Miiller elaborates. In this case 
the singular points £ £ 2 have the same ordinate. Tak- 
ing the centers of the two circles as the points £ 0 and 

to the resulting top of the profile is straight and par- 
allel with the axis of x (zero angle of attack). The 
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profiles have a certain thickness, maximum at the center, 
and the tips are not sharp-edged although they still resem- 
ble a pair of joarallel wings. Quite obviously this case 
is of no interest in the problem of wings in proximity of 
the ground. It is necessary to this end that the profiles 
have an angle of attack so that 

a = p 2 e~ 2ia 

may "be obtained. Then it may be maintained that the upper 
wing has an incidence a and the lower wing -a (wind 
in positive direction of the axis of x). The study of 
the form of the profile is, however, quite complicated. 

In any case, what interests us is the value of the 
circulation and the value of the lift; but Muller slights 
this point, confining himself to stating that the ob- 
tained lift, in contradiction to experience, is less than 
for the isolated wing. The method evidently consists in 
passing from the flow around the two circles, studied by 
Muller (reference 30), by Bonder (reference 10), and 
Lagally (reference 28), to the flow about the profile, 
and then computing the lift b^ Blasius 1 formula. It is 
quite curious, indeed, that Muller did not enlarge upon 
this method, since it is altogether simple. 

Bonder ! s findings were just the opposite (the lift 
increases to decrease with the ground distance), but the 
fact that he uses a theory which is an extension of Wit- 
oszynski 's, instead of that by Joukowski, makes it impos- 
sible to compare his results with those of the other two. 
His results, moreover, expressed in series, are quite com- 
plicated and require the use of tables compiled by the au- 
thor. Besides, Bonder ! s study has lost much of its impor- 
tance in comparison with the much more recent work by 
Tomotika. The same holds for similar studies, such as 
that by Riabouchin sky (reference 34). 

3) G-LAUERT ' § METHOD (reference 32) 

Glauert analyzes the case of the wing in presence of 
its image on an infinite plane. His method - which is 
that used by Prandtl and others - consists of taking the 
curvature of the induced flow into consideration. The 
particular feature of G-lauert*s method is the admission 
that the induced velocity must be applied with reference 
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to the midpoint of the wing, while thepoint vortex repre- 
senting the circulation around the wing must "be placed at 
the center of pressure, which is the centroid of the "bound 
vortices distributed along the chord of the wing, that is: 

c m - - e c p 

( C m = moment coefficient about the leading edge with r e~ 
spect to the center of pressure 01. The coordinate of the 
midpoint of the wing with respect to the point at which 
the point vortex is assumedly placed, results in: 

x - (0.5 -6)1 

Then with h = distance "between the two wings, a sup- 
plementary incidence is the result - the effect of which is 
additive to that of the curvature. 

Finally, the moment coefficient is written in the form 

°m = » - I °p 

which gives: 

°P - °P 0 = 4 (h) S (C P + 2 ^ (9) 

or else, reverting to the circulation (obtained "by multi- 
plying with J T 0 I) 

F - r ° = 1 (h) 2 (r + » V ° ° (9?) 

It is readily seen that this formula is almost iden- 
tical with that proposed "by the writer. Paragraph 3 of 
reference 24 reveals in fact with some simplifications: 

/ c 1 I 2 

f v a ; note that h represents the distance "between 

the two wings which in this article is expressed "by 2h; 
c f = tt ) : 

r - r o = m Os - a ) r - is* 2Cm ' v ° 1 

The equivalence follows from C m = 2c m 1 v 0 the mo- 
ment coefficient (constant) with respect to the centroid 
indicated by ~\± by G-lauert; the equivalence, however, 
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is limited to very small a. In addition, it is of inter- 
est to note that in passing from Cp to p. the change 

in the horizontal component of the velocity is completely 
disregarded; this omission, nevertheless, is not justified 
and may lead to erroneous results. In the case analyzed 
"by G-lauert, it is legitimate. 

Rosenhead, in his study (reference 18), analyzed the 
lift on a flat plate between parallel walls; hut the case 
of the off-center wing, from which he deduces the case of 
the wing with interference, seems to he afflicted with er- 
rors (reference 25a). 



C. REPORT BY I TIRO TAN I (reference 20) 



This study extends Birnbaum's method to include the 
analysis of a thin wing (two-dimensional problem) . 

Assuming the wing with a chord equal to 2, the vortic- 
ity is supposedly distributed chordwise according to the 
law : 



X 8 + C X 71 - X 1 



The self-induced orthogonal velocity becomes: 

v' = - | a - \ tee + c C| - § * 2 ) 

and the induced velocity at a point x, y of the plane 

where 

q x v x 

The coefficients a, b, and c are determined on the basis 
of 

v 1 + v 1 *v sin a = 0 
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where a = angle of attack and v y * is computed as veloc- 
ity induced by the reflected wing on the objective wing. 



we can compute the values of n in fxunction of the ratio 
2/h, where h is the distance of the wing from the plane 
(in general, we say of t/h, where t = wing chord), 
which the author represents in two charts: one giving n 
versus t/h for divers a, the other giving n versus 
a for various t/h. 

To pass from the two- to the three-dimensional prob- 
lem (wing of finite span), the author computes the self- 
induced velocity assuming elliptic lift distribution and the 
velocity induced by the free vortices of the reflected on 
the real wing on the "basis of a uniform distribution (al- 
most exactly like that adopted by the writer, reference 24)} 

The first is known to be constant; from the second, 
median along the span, follows: 



Po sing : 




dx = 



~a + ~rb = 2ttv sin a n 




Finally, put t i ng : 



1 + 



2t 
b 



m = 





a 



+ i £ 



2 v sin a 



4 v sin a 



we have 



r = Trtval mn 



where I takes account of the auto induct ion , m of the 
mutual induction, and n of the curvature of flow due to 
the induction. 
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Lastly, passing to the lift, we write: 

P = p w r (i - 



hence 



n ( 1 — — aim n j 
V 4h / 



Tani ! s method is, as already pointed out, very much 
like Pistolesi's method; the only difference consists in 
the method of computing the effect of the flow curvature* 



IV. RESEARCHES BY SUSUMU TOMOTIKA AND OTHERS 



These researches, also ci 
wyler*s report, are published 
Aeronautical Research Institut 
versity (reference 23). They 
mal transformation similar to 
in his study (reference 18): 
the flow is at first assumed t 
plane f = cp + i \|/ , shown in 
plane f the cut C G G f C r , 
sented at mid-height of t "by 
toff el tran sf ormati on : 



ted in the appendix of Dat- 
as Report No* 97, of the 
e, of the Tokyo Imperial Uni- 
are patterned on the conf or- 
that employed by Rosenhead, 
Given the plane z, provided 
o be ir rotational , and the 
figure 12, and applying in 
the plane f may be rep re- 
means of the Schwar z-Chr i s- 



— . M 
dt 



t 2 - "b s 



/(t 2 - c 2 ) (t 2 - g 2 ) 



where M is a constant. 

Introducing then the function p of Weierstrass, of 
periods 2(JUi, 2<JU 3 , where W% > 0 and (ju 3 i > 0 , the 
t ran sf ormati on 

t 2 f p ( s ) - e 3 

transforms the semiplane t in the rectangle of the plane 
s shown in the figure. 



From the preceding formula, follows: 

ds 1 

v [p(s) - e x ] [p(s) - e 2 ] 
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and 4~ ■ M Cp(s) - p(n)J 

as 

which, integrated, gives: 

f = - M [£(s) + p(n) s] 

Since the circulation is assumed to "be zero, f has the 
period 2o)i, and that leads to the relation: 

p<»0 = - Wi 

The constant M is determined so that f differs 
from i \j/ 0 in passing from s = uj 1 to s = (ju 1 + uu 3 , 

and then we have: 

f = JkM [{•<») + p(M.) s] 

Then the inside of the rectangle is transformed in an annu- 
lar zone in plane z (fig. 13) by means of the transforma- 
tion: 

s ■ 0) x + w 3 - log Z 

With w &M the velocity conjugated to v x - i Vy in 

plane 2, with w as its modulus, and with H as the 
angle which it forms with the direction of axis x, we 
have : 

fg ■ . w = | w(e - iH 

dz 

or, with Q = g + i log |wl 

df _ ^-ifi 

dx - 5 

Now the direction H of the velocity on the "boundary 
line and on the plate is known; thus, the problem in plane 
2 reduces to the determination of the analytical function 
Q(Z), provided that the real part on the two enveloping 
circles assume the prescribed values. 

Then, on the "basis of a formula established by E. 
Villat in 1912,* we have: 
* 

Le probleme de Dirichlet dans une aire annulaire, Rendi- 
conti del Circolo Matematico di Palermo, 33 (1912), 
pp. 134-175. 
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i(JU 

fi(Z) = 



iu). 

TT 



lif * (6) £,(|J io« « - » e) de 

— log z — - 6 ) d6 + iC 

ITT TT / 



r 11 /w 



where $(8) and ^(8) denote the values assumed "by f) 
on the outside and inside periphery, expressed in terms of 
angles in the center 8. The two functions 3> and com- 
ply, moreover, to the condition of uniformity: 



\|/(6 ) d6 



p 2TT 

/ 0(8) d8 = 

do do 
Applying Villat T s formula to the particular case of 
the problem, we obtain the expression df/dz and that of 
dz/ds, which for the sake of brevity, is omitted. 

The integration of the latter is quite a complicated 
problem which the author resolves "by adding the relation: 

. = - flkfll A(.) * O 0 

TT U 



where 



A(s) = - - 



s -t [a + 6 

TT 



) 



2T) 



P (! + »> 



With this function the author supplies another complex se- 
ries development. Indicating with 2a the width of the 
plate (2a =1), we have: 

2a = SJUE! {A(s 3 ) - A(s 4 )\ 

77 U ^ J 



s 3 = ttJj 



+ % - — e 3 



S 4 = tt)i + (JU 3 - ~ ^4 



TT 



8 3 and 84 being the angles at the center of the two pe- 
ripheries which correspond to points A and A*. 



The author finds atJb.er expressions for z and a# 
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The flow defined hy f from which we deduce the 

conjugated velocity s 4^ has, however, two singu- 

dz ds dz 

lar points in A and A*, where the velocity "becomes in- 
finite. Superimposing a circulatory flow of proper cir- 
culation k the infinity at A can "be removed according 
to Joukowski. With f * as the complete potential of the 
circulation flow, and with X = f + f 1 as the complete 
potential of the "built-up flow, the function of the veloc- 
ity re suit s : 



dX 
dz 



ax 

dZ; 



48 ^ = AH z --/- 

/\ds izj ~ cju x dZ/c 



dX /d z 
ds 



Since dz/ds "becomes zero in point A, dx/ &Z also "be- 
comes zero in the corresponding point, or for 

Z = e i6 4 

Now the circulation in plane Z is simply expressed "by 



hence 



*• - - S io g Z 



X = f - AJ< log z 

cTT 



Joukowski ' s equation gives: 

0 for 



df _ iK 1 _ 
dZ 2tt Z 



Z = e i6 4 



which readily gives: 



2 



4f1 



or 



4\J/ 0 u^ 2 



TT 



and in consequence: 

d X _ 2\|/ 0 u) 2 



ds 



TT 



L 



ds Js=s 4 
p(s 4 ) - p(p.) 

p ( s ) - p ( s 4 ) 



The aerodynamic action, hased on Blasius ! formula, 
results in: 
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P - i P 

x * r y 



dz 



" (S3 



C "being a closed surface containing the flat plate; or, 
after several transformations: 

P - 1 p = - * P & r mi is 

x y 2 tt / Vds/ dz dz 

'-2 

where E is a surface surrounding the circle inside of 
plane Z. The calculation, rather tedious, gives the fol- 
lowing results : 



0 



7 2 sin 3 



©3 
2tt 



3 \2tt/ 



where 3 is the angle of attack (fig* 12). 

This last is the expression for the lift of the flat 
plat e. 

The author demonstrates subsequently that when the 
distance I "between the center of the plate and the wall 
"becomes infinite, the preceding equation tends toward the 
known value: 

P 0 = 2tt a U 2 p sin 3 

Taking into consideration an expression previously found 
for a, we finally have: 




(10) 



The calculation of 6 3 and 6 4 bases upon the two 
following equations : 



e 3 + e 4 = 2p 



(ii) 
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The appended tables give the values of 8 ; 



function of g and 
we may then compute: 



q = e 



o>1 nl 



2a 
H 



2 

sin 8 



V2tt 



) 



St) 

2tt/ 



. Knowing 8 3 

3 -,2 
- 6 

2 



an d 8 4 in 
and 8 4 , 



t3 



3 

2-nJ 



(13) 



Formulas (10), (11), (12), and (13) are fundamental for 
the calculation. 

The author deduces from (12) the following series de- 
velopment : 

cos \ (6 3 - 6 4 ) = cos p [- 2q + 4q 3 - 4q 5 (3+2 cos 2 p) 

+ 32 q 7 (1 + 3 cos 2 8) 

- 2q 9 (37 + 324 cos 2 p + 80 cos 4 P) 

+ ] 

This formula combined with (ll) is the "basis of his append- 
ed tabulations. 

The values of P/P 0 and of 2a/H (or l/H) deduced 

from (10) and (13) are included in the tables and plotted 
in figure I4« 

Thus it is readily seen that for small 2a/H, the 
ratio P/P 0 is always less than 1, regardless of angle of 

attack; but for small 3, the values of P/P 0 become 

greater than 1 when 2a/H exceeds a certain critical value 
which is an incremental function of p. For medium an- 
gles of attack the lift is manifestly increasing in accord 
with experience, by small distances. 

Tomotika f s calculations do not include very small dis- 
tances Ht In particular, they do not come to an end at 
the cancellation of the slot s between wing trailing edge 
and the plane of infinite length. Obviously, this condi- 
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2a 2 

tion prevails when H = a sin P, or for — - = — - — 

H sin p 

In ^articular, % = 6.47 for p = 18°, and ™ = 

0 P 

3.4 for p = 36 . Compared with Datwyler ! s data, the 
results seem to fairly agree, hut not enough to remove 
all doubts. (Fig* 16 indicates the value of P/P 0 ob- 
tained "by Datwyler by a small circle; that is, 0.630* ) 

Lastly, Tomotika gives an approximate development 
of P/P 0 in powers of 2a/H, and compares the results 
with that of the rigorous calculation (applicable to small 
values of 2a/H). The expansion for P/P 0 in powers of 
2a/H finally gives: 

P_ = i - Jin_g fg^ + 1 (4 _ 3 cos 2 p) (gg) m 
P 0 2 V H / 16 v v H / 



4 

+ (32 - 57 cos 2 P + 22 cos 4 P) 

512 V H , 



The lift calculation, evolved on Cisotti's ideas, is 
given in appendix I. 

In another report (reference 25) Tomotika attacks the 
problem of lift on a flat plate in a stream bounded by an 
infinite plane wall. 

Although these problems do not seem to have any prac- 
tical significance, they are, nevertheless, of interest 
from the theoretical point of view. The result is shown 
in figure 15, from which it is readily seen that there is 
an increase in lift, increasing with the angle of attack 
and ratio 2a/S (decrease with distance from the plate). 

The case in which the interval s between leading 
edge and plate is zero, in Datwyler r s report, leads to a 
negative lift, as is readily understood when observing 
that in his treatment there is always a force toward the 
wing. Datwyler*s problem may be considered the limit of 
Tomotika 1 s problem. 
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V. APPROXIMATE FORMULAS DERIVED BY PISTOLESI 'S METHOD 



It might he of interest at the conclusion of this re- 
view to check whether my own method suggested for the cal- 
culation of the infinite hiplane lends itself to conven- 
ient application to assure a high degree of approximation 
in the solution of the prohlem of ground interference. 

The concept of this method may he hriefly summed up 
as follows; 

a) Compute the circulation ahout the wing on the "ba- 
sis of the velocity (further induced asymptotically) pos- 
sessed "by the wing at a point l/4-chord length from the 
trailing edge (point K ) ; the lift is obtained by multi- 
plying by c f I the component of this orthogonal velocity 

1 d C-v-) 

of the axis of the profile (c f = -r t — ; for airfoils with 

s 2 da 

small camber, c ! = tt ) . 

b) Substitute, for the effects of induction and lift 
calculation, the wing with a vortex located in the center 
of gravity of the circulation; or else in a point situated 
I I length from the leading edge (point I), with a proper 
doublet added in the same point. In the case of a flat 
wing, the doublet is zero and the wing may be replaced by 
a simple vortex located at the said point. 

Now let us see how this concept can be applied to the 
wing in the presence of its reflection; we limit ourselves 
to the case of a flat wing. The more simple application 
is that disregarded in the induction calculation: the an- 
gle of attack of the wing (fig* 16). 

YIe have then at point K: 

| 

r 2 

v n = v 0 sin a + — — 



and consequently: 



2TT 4h 2 ♦ £ 



i 3 



r = tt i v„ = r n + r 



n ~~ 0 2 7 2 

4h + — 

4 



which givcsJ 
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r " r ° C 1 + 158*3 <14) 

The velocity in correspondence with point I is 

V ° " 4tt h 



and finally, the lift is: 

' = ''('•■ its) 



or 



P = p tt I v 2 sin 



a 



i - - l - (i + 

4h V, 



1 ^ 

lih 5 "; 



sin a 



1 + 



Po 



1 - 



4h 



(l + At ) sin a ] [i + JL 
> 16h > J L 161i J 



Applying this formula to the case of a = 4^°, we find: 

0.1 0.2 0.3 0.4 0.5 



(15) 



JL = o 



4h 



1- = 1 



1.002 1.023 



1.062 



1.118 



1.188 



Plotting these results (fig. 14, dotted line), they 
are seen to remain consistently above the rigorous values 
of Tomotilca, and not widely at variance at high incidence. 
As regards the inclination of the wing in the calculation 
of the induced velocity, the approximate formulas (3) and 
(8) of my previous report could he applied, "but other for- 
mulas more rigorous and no more complicated than they, are 
preferable. We have (fig. 17): 

r 

v n = v 0 sin a + sin Y 

27Td 

p 

r = tt I sina.+ sin Y) 

0 2-nd 

= tt 1 v A sin a + Ulu (sin Y_ cos a - cos Y sin a) 
0 2tt d 0 



Then 



„ f „ v _ I co s a, 
sxn Y n = 

0 2d 



v 2h 
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so that 

r = r + --2- f~ cos 2 a - 2h sin 

o 2d \2 J 

or with 

after simple changes: 



d 2 = 4h 2 + \- cos 2 a 
4 



1 + 



(k1 



cos a 



b£M 

1 + 4h Sin a 



(16) 



where k is the multiplication factor of P 0 • We have 
also: 



P = p r 



v 0 - 



2h + ~ sin 



) 



which, after various changes* gives: 

P 



with 

whereby, again: 



p 0 k» 



k' = k (1 - kr) 



(17) 
(18) 



k = - 



cos a 



(19) 



1 + -fr sin a 
4h 



*In taking account of the wing curvature ( c m x ± 0), we 
arrive by similar treatment at a result which is always 
exiDr es sibl e "by (18), with 

c 

2tt 

sin a (1 + X sin a) 



k = 



sin a ( 1 + X 2 cosa) - -p~ ^ 2 



X sin a 



r = TL 



1 + X sin a 



(190 
(20 r) 



whereby A = 



4h 
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— — sin a 

t -_ (20) 

1 + -T7- sin a 
4h 



Applying this formula to the case of a = 4.5°, we obtain 
the following tabulation: 

JL m 0 0.1 0.2 0. 3 0.4 0.5 

4h 



JL « 1 0.9943 1.0085 1.0404 1.086 1.146 
?o 

These figures are equally superior to those of Tomo- 
tika, hut much more approximate than the former, especially 
for small l/h (up to 0.8). They are shown in figure 14 
as dot-dashed line. 



A markedly improved approximation is obtained if the 
calculation of the induced velocity on the basis of a sin- 
gle vortex concentrated at I* is abandoned in favor of 
the field of the velocity around the reflected wing, for 
which we proceed as follows: 

The velocity about a flat plate of length 2, assumed 
with a velocity of the components v x and v 2 and with a 
circulation T is: 



i - - 



- 1 



2tt y~ 2 a - i 



The application of the above formula to the field 
about the reflected wing, gives: 



- v. 



cos a 



v 2 = - v c 



sin a 



and w e have : 

w = v 0 cos a + i v 0 sin a ~p~=~ ~ 1 J~ "7=^= ^ 21 ^ 

V z -1 ^ fT v z -1 

For the determination of the circulation, it is nec- 
essary to determine w at point K on the real wing 
which, relative to the reflected wing, has the ordinates: 

x K a 2h sin a + § cos 2a 
yg e 2h cos a ~ § sin 2a 
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In general, denoting with X the ratio of distance 
h and wing chord I , X = h/l , the previous formulas can 
he written as: 



y = 4x cos a - 2 cos 2a 



(22) 



To determine the lift - whatever the circulation - 
we must compute w at point I of the real wing which, 
relative to the reflected wing, has the ordinates: 



Xj = 4x sin a - J cos 2a 

yj = 4x cos a + § sin 2a 

From (21) follows: 

v ■ v cos a - v sin a (Bx + Ay) + — 1 
x o o 2tt 

r 

v y s - v 0 sin a (Ax - By ) + ™ A 



(23) 



(24) 



where 
A = ± 



B = ± 



f\ (x 2 -y 2 -l) + \ J ( x 2 - y a-l) r +""4x2 y« 
/- |(x 2 -y 2 ~D + i/(x 2 -y 2 -l) 2 + 4x 2 y 2 

77?-7^7"~4x 2 ""; 2 



(25) 



With regard to the signs of the radicals, it is readily 
seen that A will have the sign of x; B, the sign of ~y. 



The normal component of the real wing at K will he 
.n 2a + • 

nates of K: 



v x sin 2a + v y cos 2a, where x and y are the coordi- 



v n =v Q |[cps a-sin a(Bx+Ay)] sin 2a-sin a(Ax-By) cos 2a) + 



+ — (B sin 2a + A cos 2a) 
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which, since F = 2tt v n gives: 

p = 2tt v [cos q-sin q(Bx+Ay)] sin 2q-sin q(Ax-By) cos 2a 
0 1-B sin 2q-A cos 2q 



or 

T 1 - ( BgXg+A^-y ) sin 2q - ( ■ A -K x K"* B K y K ~ !) cos 2q 
T 0 1 - Be sin 2q - A& cos 2q 



(26) 



To pass to the value of lift P, it is necessary to 
know the hori zontal velocity corresponding to point I, 
which will "be v x cos q - Vy sin a: 



=v Q | [cos a-sin q(Bx+Ay)] cos q+sin q(Ax-By) sin qj> + 



p 

+ — (B cos a- A sin a) 
2tt 

or, writing T = 2tt v Q k sin a (where k ■ r/r o ) and 
making several changes: 



IftI 
v o 



1 - sin q <f[ Bj (xj-k) + Ajyj] cos a - [Aj(xj-k) + 

+ Bjyj - l] sin qj (27) 

which, for ~ ss r, finally gives* 
v o 

P = P 0 kr (28) 

The calculation is quite rapid since the quantities 
A and B are easily computed, after the expressions x 2 •* 

y 2 ~ 1 and 4x 2 y 2 have "been calculated. 

Applying the method to the Case of a = 18° (i.e., 
very high incidence), and for the whole range of values 
X (or of X a l/x), analyzed by Tomotika (X = 0 to 3)), 
the agreement with his test results la practically perfect, 
so that the two curves coincide. This is valid for the 
plane problem - that is, the problem of the infinite wing. 

To analyze the problem of the wing with finite span 
necessitates taking into account the velocity induced by 
the free vortices and the calculation is notably compli- 
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cated. In view of this complication, I felt justified 
in my study (reference 24) to confine myself to consid- 
ering this effect "by the method of horseshoe (or trail- 
ing) vortices (lift distribution by uniform law); it is 
probably possible to improve the research method "by 
adopting Ferrari f s method (reference 36) for the calcu- 
lation of the biplane of finite span, which is based upon 
elliptical lift distribution. Still, even supposing it 
is possible to substitute the wing with a distribution of 
vortices on the aerodynamic center of the wing and with a 
distribution of doublets on the central axis of the pro- 
file, its applicability with good approximation is lim- 
ited to moderate ratios of l/h. 



Translation by J. Vanier, 
National Advisory Committee 
for Aeronautics. 
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Fig. 1 




x--. Biplane alone 

o Plate 38 mm below lower wing 

™{~ Reflection 76 mm between lower wings 

Airspeed 12.19 m/sec 

Wings of biplane 76 x 457 mm 2 



Figure 1. 
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Figure 5.- Lift of flat plate in 

proximity of the ground 
(partly approximated). 
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Figs. 3,4,8,9 
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Figure 3.- t plane 

(original plane). 
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Figure 4.- Lift determinar 
tion. 



2s = 5 10 20 mm 




0 4 8 0 4 8 

a ,deg. a ,deg. 

Figure 8.- Zurich experiments Figure 9.- Comparison "between 

(reflection method) theory and experiment. 
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Figs. 12,13,14 
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, Z plane 
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H Figure 13.- Tomotika's 

!! experiments 
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i Plane Figure 14.- Tomotika's results 

Figure 12.- Tomotika's experiments. compared with those of 

Pistolesi and Datwyler. 
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Figs. 15,16,17 




Figure 15.- Toraotika's results 

for the case of 
the plane above the wing. 
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